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Abstract. When particles are multiply scattered by a random poten- 
tial, their momentum distribution becomes isotropic on average. We 
study this quantum dynamics numerically and with a master equation. 
We show how to measure the elastic scattering time as well as charac- 
teristic isotropisation times, which permit to reconstruct the scattering 
phase function, even in rather strong disorder. 

1 Introduction 

These days, ultra-cold atoms permit to study disorder physics in a quantitative man- 
ner. In particular, quantum transport of non-interacting particles can be studied in ID 
PQ, 2D |2|3|4j and 3D |5|6j . One of the latest development, suggested in [7], consists in 
launching a quasi-monochromatic wave packet inside the bulk of a random potential 
and monitoring the single-particle momentum distribution |3|4| . The long-time evo- 
lution leads to coherent backscattering (CBS) and coherent forward scattering (CFS) 
signals, which are linked to weak and strong localization, respectively |5]. In this pa- 
per, we analytically and numerically analyze the short-time dynamics, which contains 
valuable information about disorder strength and scattering characteristics. In par- 
ticular, we show how one can not only measure the scattering and transport times, 
but also reconstruct, by an angular Fourier analysis of the momentum distribution, 
the complete scattering phase function. 



2 Momentum-distribution dynamics 

Starting from the momentum distribution rifc(O), we seek to compute rik(t) = |^fc(^)| 2 
at later times t under the evolution with H = H^ n + V(r). The over-bar {.} denotes 
an ensemble average over realizations of the random potential V(r). 

We simulate this evolution numerically by solving the two dimensional (2D) Schro- 
dinger equation for matter waves with H^ n = p 2 /2m using a finite- difference method. 
The initial wave function has a small width hAk around a finite average momentum 
hko- The random potential V(r) is an isotropic speckle with a Gaussian spatial co- 
variance length a and amplitude V. The natural energy and time scales of the system 
are E a = h 2 /ma 2 and t a = ma 2 /h. Figure [l] shows numerically calculated momen- 
tum distributions, averaged over 1000 disorder realizations, for two different regimes: 
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Fig. 1. Density plots of momentum distributions from the numerical simulation. Upper row: 
Experimental regime of g], where ko = 0.9/cr, V = 0.35^, and Ak = 0.07/co- Bottom row: 
ko = 1.8/ a, V — 0.7 E a and Ak — 0.035/co- Parameters are chosen such that the weak- 
disorder ratio A — 2V 2 / (koa) 2 « 0.3 of perturbation theory is the same 9 . Right panels: 
the peak value rik of the initial momentum component decreases exponentially, Eq. |5]). 
The indicated fits yield the elastic scattering times Jq 1 = r s = 6.7t a and 2At a , respectively. 



rather isotropic scattering with koa = 0.9 (upper row) and rather anisotropic scat- 
tering k^a = 1.8 (bottom row). 

As the initial peak at +ko decays, due to elastic scattering by the random po- 
tential, other momenta are populated on the (disorder-broadened) Rayleigh ring of 
radius fco- Eventually, a CBS peak appears at — fco- In the remainder of this paper, 
we rather focus on the dynamics that produces the isotropic background. 

From an analytical point of view, since unitary evolution conserves energy, it is 
appropriate to study the joint distribution rik(E,t) of energy and momentum, with 
marginal rik(t) = J T^rik(E,t). The density propagates according to 



n k (E,t) 



duo 
2^ 



(1) 



where the intensity propagator (E,uj) obeys the integral equation of motion |9|10j 



■-y k {E)]$ hk ,(E,u>) = A k (E) 



Jkk' 



^U kp (E,uj)$ pk >(E,uj) 



(2) 



This form is valid for small frequencies uj <C E and can thus describe the dynamics 
at long enough times t ^> h/E. The spectral function A k (E) = 2tt (k\ S(E — H) \k) 
is the probability density of a plane- wave state k to have energy E. The vertex 
U kp (E,u) describes all scattering events that couple i/j and It is related to the 
elastic scattering rate j k (E) by the so-called Ward identity |9|llj 



7fc(£) = ^A p (£)/7 pfc (£,0). 
p 



(3) 
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This relation is a variant of the optical theorem, and expresses particle-number con- 
servation. Using (—iuj)<P kk '(E,uj) from ([2| in ([!]) then yields the equation of motion 

d t n k {E,t) =A k (E)n k (0)S(t) - lk {E)n k {E,t) 

+ A k (E) J ^e-^^C7 fcp (^,a;)^(E,a;)n fc ,(0). (4) 

Initially, the populations at k' ^ ko outside the initial wave packet are small. 
The most noteworthy evolution is the decay of n ko as described by the first line of 
@, whose solution immediately leads to n ko (t) = 0(t) [J ^A ko (E)e~^ h o^ t ] n ko (0). 
This is a superposition of simple exponential decays, weighted by the spectral density. 
As function of E, A k (E) is rather sharply peaked around the effective excitation 
energy E k , whereas / y k (E) is smooth. Therefore, the dominant behaviour is 

n ko (t) =e-T° t n fco (0) (5) 

with the elastic scattering rate 70 = 7fc (E , fc ). 

The corresponding numerical data are shown in the right-hand panels of Fig. [T] 
For the two regimes of rather strong disorder studied here, we find r s = 6.5t a for the 
first case and r s = 2At a for the second. Lowest-order perturbation theory |9|12j over- 
estimates the scattering strength and predicts r s = 2M a and r s = 1.4t a , respectively. 
We have, however, verified that if the simulation is run for V 2 /EE a = 0.005 <C 1, 
the measured value of r s agrees with |9|12j . 



3 Diffusive isotropisation 

3.1 Master equation 

Equation Q describes the equilibration of populations due to potential scattering. 
Whenever the potential-scattering vertex U kk >(E,uj) is regular for uo — ^ 0, the singu- 
lar ^-dependence of <& kk > (E, uj) dominates the cj-integralQ Using ([!]) and the Ward 
identity one then arrives at a Markovian master equation of the Pauli type: 

d t n k (E,t) = J2 Uk P^ iME)n p (E,t) - A p (E)n k (E,t)} . (6) 
p 

Assuming a factorised solution n k (E,t) = A k (E)n k (t), we then obtain a closed mas- 
ter equation for the momentum distribution, d t n k (t) = ^2 U kp [n p (t) — n k (t)], with 

transition rates U kp = J T^A k (E)A p (E)U kp (E). The spectral functions fix the ener- 
gies at E p w E k . Consequently, we can focus on the angular probability density n(0) 
on the circle k — /co(cos#, sin#) of elastic scattering. There, scattering from 6 to 0' 
is described by the on-shell value U{6 — 6 r ). It is helpful to re-introduce the elastic 
scattering rate 70 = ^0 / d0U(0), in accordance with ^0 = v(E ko ) denotes the 

1 U(uo) w U(0) can only be factorised from the cj-integral in Q if it has no singular 
dependence in uj. This is true at short times, since the single-scattering term does not 
depend on frequency at all and the frequency dependence is smooth for the first interference 
corrections. But the CBS contribution develops a diffusion pole at uj = and backscattering. 
Thus, the momentum distribution retains a memory of the initial condition and displays the 
CBS peak [7 . Even more singular terms appear at the onset of Anderson localisation and 
produce the CFS signal [8]. These phenomena are beyond the master-equation approach 
developed here. 
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Fig. 2. (a) Fourier coefficients of the momentum distribution on the elastic scattering circle, 
up to m — 4, as function of time, for the case koa = 1.8. Straight lines are fits to the 
exponential decay, Eq. for t a < t < bt a . (b) and (c) Fourier coefficients u m of the 
scattering phase function, as measured from the numerical data (full red/blue), and at the 
single-scattering approximation, Eq. ^ (dashed black). The insets show the reconstructed 
phase functions, using moments up to m = 4. 

density of states. Finally, the elastic isotropisation of momentum is described by [see 
also Eq. (6.83) of [13]] 

p2tt 

dtn(0, t) = 70 / d0'u(0 - 0') [n(0', t) - n(0, t)) , (7) 
J o 

where u(f3) = U((3)[f dOU(6)]~ 1 is the normalised differential cross section of elastic 
scattering, the so-called phase function. 

3.2 Characteristic times 

The convolution Q on the unit circle can be readily solved by Fourier analysis: the 
coefficients c m = jd0e %rne n(0) obey the decoupled equations of motion 

d t Cm(t) -70 (1 - U m )c m (t) = -T^Cm(t) (8) 

in terms of the phase-function Fourier components u m = J dj3e Lm ^u(j3). By virtue 
of uq = 1, the zeroth mode cq (i.e., the number of particles) is conserved. The first 
harmonic c± decays exponentially on the transport-time scale t\ = r s /[l — (cos/3)], 
with the usual notation (cos/3) = u\ = J d(3 cos f3u(j3) [9 . Higher harmonics decay 
with their own times r m = r s /[l — (cosra/3)]. 

Only for a completely isotropic phase function with (cosm/3) = all times r m 
are strictly identical to the elastic scattering time r s . Things are more interest- 
ing for spatially correlated potentials. Indeed, at the single-scattering approxima- 
tion, the Gaussian correlator uj^J ex exp{— a 2 (k — p) 2 /2} has the phase function 

u^\/3) = exp{k 2 a 2 cos f3}/[27rIo(k 2 cr 2 )] where Iq(-) is the modified Bessel function. 
The corresponding Fourier coefficients 

U W = I m (k 2 a 2 )/I (k 2 a 2 ) (9) 

decrease as a function of the order m, such that the decay rates r" 1 increase. In 
other words, the highest Fourier modes disappear first, giving way to an ever-broader 
background. 

In Fig. [2(a), we show the time evolution of the first four Fourier coefficients c m 
of the angular distribution in the case ko<j = 1.8, as obtained from the numerical 
simulation data shown in Fig. [l] by a radial integration over a small interval around 
|fe|=*o. 



Will be inserted by the editor 



5 



3.3 Reconstruction of the phase function 

From the measured decay of c m (t), one can compute the Fourier coefficients u m via 
Eq. (J8J), and thus reconstruct the phase function u(/3). Fig.[2|b) and (c) show the u m as 
measured from the numerical data, together with the single-scattering approximation 
([9]). In the first case ko<j = 0.9, the higher-order Fourier coefficients of the momentum 
distribution are very small. Consequently, the scattering phase function, shown in 
the inset of Fig. [2jb) , is rather isotropic, as expected for k^cr < 1. In the second 
case koa = 1.8, the scattering is much more pronounced in the forward direction, as 
expected for ko<j ^> 1. 

Compared to the bare single-scattering phase functions (in dashed), the full re- 
constructed phase functions show a tendency towards enhanced backscattering. We 
may attribute this to weak localization corrections associated with short scattering 
paths, which are only known for isotropic scatterers [T4] . 

4 Conclusion 

We have analyzed the momentum isotropisation of a quasi-monochromatic wave 
packet inside a spatially correlated random potential. Results from the numerical 
calculation are well reproduced by a master-equation approach, even in regimes of 
strong disorder. We show how to reconstruct the phase function of elastic scattering 
from the measured exponential decay times of the angular components of the mo- 
mentum distribution. This method is directly applicable in present-day experiments 
[314] . and can serve as an in situ calibration of random potentials, for instance when 
investigating the impact of disorder in the presence of interactions |15|16j . 
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